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Abstract 


A nonlinear wave mechanical equation is proposed by inserting an 
imaginary quantum potential into the Schrödinger equation. An explicit 
expression for its solution is given under certain assumptions and it is 
shown that it entails attenuation related effects as non-unitary evolution, 
non-exponential quantum decay and entropy production. In the quantum 
hydrodynamical formulation the existence of circulation effects for the os- 
motic velocity field is established. Finally, a time-invariant equation for 
the probability density is derived, analogous to the tensor Lighthill equa- 
tion in aeroacoustics, which admits both retarded and advanced solutions. 


1 Introduction 


One of the most spectacular aspects of quantum theory is that it can be formu- 
lated in a variety of ways that even appear not to be equivalent at first sight. 
In some occasions the mathematical formulation is identical but the conceptual 
framework differs significantly, in other instances even the mathematical formal- 
ism varies. It is well-known that soon after the inception of quantum mechanics 
de Broglie proposed an alternate formulation based on a guiding wave field, 
an idea rediscovered later by Bohm [I]. Almost at the same time it became 
possible to formulate the new theory in terms of a hydrodynamical model [2] 
and this line of enquiry was later continued and perfected by Takabayasi [3]. In 
both formulations, although in a different context, a prominent place is given 
to what is called quantum potential. The Bohm quantum potential should not 
be considered as another classical potential function acting on physical systems 
[4]. It is a rich physical concept with broader mereological implications and dis- 
plays with clarity the inherent wholeness and non-separability of an interacting 
many-body system [5] since it does not necessarily diminish with distance and 
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it is manifestly non-local. On the other hand its very reliance on a probabilistic 
measure, namely the curvature of the amplitude of the wavefunction [4], seems 
problematic within the scope of a deterministic theory like the de Broglie-Bohm 
one. In the next section we examine some fundamental considerations, introduc- 
ing a non-linear wave mechanical equation and relate its solution with the linear 
Schrodinger equation. In section 3 we examine the hydrodynamic formulation 
and examine the effect of the quantum potential to the osmotic velocity. A time 
symmetrical wave equation with a non-homogeneous term for the probability 
density is also derived that applies to the unmodified dynamical problem. 


2 Foundational issues 


Here we modify the Schrodinger dynamics by introducing an imaginary Bohmian 
potential and then make a polar substitution for the wavefunction. In that way 
we have introduced a quantum potential term in both the Hamilton Jacobi and 
the modified continuity equation as a point of contact and consequently we 
recombine those two expressions in a different manner assuming at the same 
time that the phase of the nonlinear wavefunction is approximately the same 
with that of the linear equation. Our point of departure is a modified nonlinear 
wave equation of the form 
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the Bohm quantum potential discussed earlier. Making a standard polar substi- 
tution of the form ® = Re’S/" and S the multivalued phase of the wavefunction 
gives the continuity equation 
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which can also be written as 
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The Hamilton-Jacobi equation reads as 
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It follows from (4) that even if the phase does not exhibit any spatial variation 
the probability density may still be time-dependent. It can be reduced to a 
nonlinear equation, transformed to a diffusion equation like (3) and be readily 


solved provided that £ is positive, otherwise the coefficient of diffusion will be 
negative and the problem ill-posed. On the contrary € < 0 describes amplifi- 
cation effects. A diffusion equation requires initial conditions to be imposed in 
order to deal with a well-posed problem, in other words if we use final condi- 
tions for (3) or (4) this engenders the possibility of a blow-up to infinity and a 
singularity in a finite time interval. This means that there is a preferred time 
direction and it is not possible to retrodict solutions in the backward time di- 
rection. If we allow only the first order derivative phase term to remain we 
obtain a diffusion equation with convection that can still be solved exactly. The 
existence of the Bohmian potential in (1) allows us to establish a connection 
between the Hamilton-Jacobi and continuity equations as 
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The right hand of the above expression, which can been seen as a source term for 
the continuity equation, is related to the classical Schrödinger equation [6]. The 


time-independent version of (6) is easily obtained through the relations oe = 0, 
25 = —E. By introducing the convective derivative which is expressed as 
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Set We v Oy 7 
da Ram” m 


it is possible to solve (6) for p = R? 
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As in [7| we can write 
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where L the system Lagrangian. Multiplying the last two equations we obtain 
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In that way we have generalized the expression for the quantum hydrodynamic 
propagator given in the de Broglie-Bohm theory [7]. Notice two alternative 
gauge-independent expressions are given. These are formulas relating the wave- 
function ® at two different times and spatial locations, through an exponential 
factor that involves integration over past times. It would be interesting to find 
a relation between ® and W at the same time and location. In order to do this 
we need to make an assumption about the phase. We assume that the phase 
of the unperturbed system does not change significantly when we include the 
non-linear term. It is reasonable to assume that the imaginary term will more 
appreciably affect the evolution equation of the probability density and to a 
lesser extent the time development of the phase since the Hamilton-Jacobi that 


governs the evolution of the phase equation will remain unaltered. Certainly, we 
can add correction terms to this initial phase and use (5) to obtain iteratively 
better approximate expressions. Using (8) and setting € = 0 we find 
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where Y = Ae‘S/”, Multiplying both sides with e*%/” and assuming that 
R(ro, to) = A(ro, to) we find 
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Hence we have found a proportionality relation between ® and Y. Notice that 
the quantum potential Qo corresponds to the state Y. This approximate expres- 
sion indicates that the total probability amplitude is rescaled by an exponential 
integral factor with local phase dependency. When the finite duration of the 
time aperture becomes instantaneous so tg — t the exponential becomes unity 
and we recover the linear solution V at equal times. Obviously we recover 
the linear limiting case too when ¢ > 0. Within a many particle description 
the phase modifies a coincidence probability density and features the non-local 
character of the generated dissipation. It must be remembered that an essen- 
tial hypothesis, apart from the inclusion of the imaginary Bohmian potential, 
was the assumption that the phase dependence of the modified wavefunction 
is identical to the linear Schrödinger equation. Only then can we obtain this 
compact relationship between the two wavefunctions which would otherwise be 
unrelated. Since ® is a solution of a nonlinear equation (1) in contradistinction 
to W it does not satisfy the superposition principle due to the existence of the 
scale factor. Nevertheless, it has the same phase so it may exhibit interfer- 
ence and correlation effects exactly as a coherent superposition of Schrödinger 
equation solutions. A coherent state for example will display the same phase 
properties but an amplitude with a different spatio-temporal modulation. The 
presence of the state associated quantum phase illustrates the kind of indivis- 
ible, essential unity which is characteristic of the whole quantic system, more 
clearly pronounced in the multi-particle case. Normalization concerns dictate 
that ® goes to zero at infinity and it may even involve non-Markovian evolution 
and entail memory effects [8]. Using (12) we can re-express (1) as a non-linear 
Schrödinger equation with an added integral term. This is not a particular 
useful form but it is included here for the sake of completeness 
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Effectively the existence of the imaginary Bohm potential means that the quan- 
tum action becomes complex 
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where a vector potential has been introduced in the usual manner as v = xs — 


+A. The above relation suggests the following form for the non-relativistic 
propagator from place ro, to to r,t 
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This expression clearly implies non-unitary evolution. Again it must be empha- 
sized that it is this phase or action of the system that generates dissipation, so 
it is a kind of self-induced decay and unlike the influence functional or similar 
approaches we do not need to have recourse to a heat bath that affects the sub- 
dynamics of the system under examination. This is of great consequence and 
illustrates the sharp conceptual difference between the way irreversibility ap- 
pears in the present formalism and more common reservoir-system approaches. 
In fact, it may be argued that strictly speaking it is not possible to trace out 
environmental degrees of freedom and divide rigorously a system from its sur- 
roundings, there is no sharp distinction between measured system and environ- 
ment since the quantum phase cannot be appropriated to individual particles 
but points to the irreducible quantum compound and reveals its substantial 
unicity. If the action in the frequency domain is analytic in the upper half 
then its real and imaginary parts are not independent but they are related by 
a Hilbert integral transform 
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the frequency domain representation of the integral scale factor and W(r,w) 
that of W(r,t). Furthermore using the Fock-Krylov formula [9] for the survival 
probability of an unstable system we find 
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which in general deviates from exponential decay law even if the density of 
energy states is a Lorentzian function . The decay rate is 


_ 2 gs + =-(VS — gA)? + V(r, 7) 


Or ams L 19 
h o ¥ Ja [BE trk (SGA)? +V (r,r)]ár a 


It is straightforward to generalize (1) for a many-particle system (of equal masses 
for convenience) or even a continuous one like a boson field. It is evident that 
if it is possible to decompose the phase as a sum of functions S = S1 + S2 +.. 


then the multiparticle wavefunction is separable. The rate of change of the total 
probability P = f pdx 
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We see then that the total probability is not conserved but diminishes with time. 
What is important is that (1) entails the possibility of entropy production. Using 
the entropy formula S = — f plnpd?x we obtain 
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which is positive. The positivity of entropy production is a general result fol- 
lowing directly from the modified dynamics and does not rely on the phase 
approximation made earlier. 


3  Retrocausal solutions and osmotic vorticity 


It was stated in the introduction that quantum mechanics can be expressed 
in terms of quantum hydrodynamic variables, namely the probability density 
and the momentum field, even though the two formulations are not equivalent 
as certain hydrodynamical solutions may not satisfy the Schrodinger equation. 
The inhomogeneous continuity equation for the probability density is analogous 
to the conservation of mass equation in fluid mechanics. It is possible to derive 
an equation analogous to the momentum Navier equation [8] [6] [7]. Employing 
(4) we find the continuity equation for the flux momentum 
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where the quantum force is the sum of the classical and quantum potentials 
dv 


It is possible to write (22) in a tensor form as in [6] 
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where the probability quantum stress tensor is written as 


Qui, (24) 
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We recover the familiar expression for the tensor equation [6] in the limiting 
case € + 0. Following a development of the Navier-Stokes equations found in 
the aerodynamic sound literature {Ti} we take the time derivative of the 


continuity equation for the probability density (4) and the divergence of (24) 
and by subtracting them we eliminate the convection term. In that way we 
combine the two in a new exact form as 
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which contains both second and first order time derivatives and in order to find 
its time evolution we need to know initial and final conditions for p and p 
both. The imposition of two-time boundary conditions is already familiar but 
an important remark that needs to be made is that here the primary object of 
interest is not the quantum state but instead the probability function p which 
we believe lays on a higher level of ontological density. The second term in 
(26) contains non-constant coefficients so the equation cannot be solved exactly 
but notice that it resembles the telegrapher’s equation which admits solutions 
expressed by Heaviside step functions. This indicates at least qualitatively that 
the solutions of (26) will propagate not instantaneously but with finite speed. As 
in the aeroacoustics case we subtract from both sides a second spatial derivative 
term multiplied by an arbitrary reference velocity co that states the propagation 
of p and obtain 
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where the driving Lighthill tensor which includes a quantum potential compo- 
nent with highly non-classical properties is expressed as 


(27) 
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For a two-dimensional problem we could write (27) in terms of a stream function 
for the flow velocity v. If we omit the quantum potential term and set £ = 0 in 
(27) we obtain an evolution equation invariant under time reversal which admits 
solutions in both time directions. In contrast, the Schrodinger equation admits 
solutions only in one direction in time and its complex conjugate in the opposite 
one. In that limit (27) is reduced to an inhomogeneous wave equation with a 
Lighthill source term [Ii] [12] and its solution is 
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where r = |x — y| the distance between a reception point x and the source y. 
Even though this is a formal solution it still allow us to extract some signifi- 
cant conclusions by including both retarded and advanced terms. The existence 
of the retarded term is rather uncontroversial but this cannot be said for the 
advanced term since unavoidably introduces an element of finality and retro- 
causation [13]. A similar issue is well-known in classical radiation theory where 


the advanced terms are excluded on physical grounds, but it is not clear if the 
same kind of argumentation could be employed here. Notice that in the quan- 
tum version of the single particle Lighthill equation derived here the quantum 
potential is embedded in the right hand side tensor. This implies nonlocal ef- 
fects and strongly suggests that we cannot discard advanced solutions without 
loss of essential physics. Furthermore, for two or more particles the coincidence 
probability density existing in the 3N dimensional configuration space is con- 
nected to the total probability that appears in Bell type theorems which is a 
clear indication that the advanced term needs to be included in order to violate 
the classical correlations. This is obviously not a mathematical argument and 
the matter needs further consideration but it is plausible to assume that the 
effect of advanced terms in a many-particle context would definitely affect cor- 
relations and violate the usual assumed conditions for Bell-type theorems {14}. 
Obviously in the context of a classical electromagnetic or acoustic field there 
is no possibility of violating that family of inequalities but this is cerainly not 
the case for a quantum coincidence probability. One of course could escape the 
conclusion by modifying the ordinary Schrödinger dynamics. Introducing now 
the stochastic osmotic velocity [7] [15], which is also implied in earlier formulas 
as in (21) and (25), 





nh (VO Vor 
ue (SiS) en) 
gives 
h de f* 
u = —Vin|v/? — = | (> + IWV vt vv) dr. (31) 
2m M Jig \ OT 


This is a measurable, phase dependent diffusive velocity guided by the deter- 
ministic flow velocity field. The many particle expression exhibits dependency 
on the remote site locations which is a manifestation of nonlocality. Notice that 
integration over past times can be extended in principle to the initial proba- 
bility distribution which makes a common cause explanation [16] implausible. 
The total velocity is the stochastic one in addition to the standard flux veloc- 
ity. The intrinsic statistical natural of the osmotic velocity field does not allow 
deterministic trajectories and the rule of non-crossing is not applicable. The 
dilatation of the osmotic field is given by 
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More interesting is the vorticity of this stochastic velocity field. Taking the 
curl of (31) we confirm that the first term is zero as the gradient of a scalar 
function but the integral term has non-zero contribution and some interesting 

implications. Using the fact that V x v = 0 yields 
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We see then that the quantum velocity has a non-zero component related to 
circulation which could affect in principle the electron current in a conducting 
ring as in the case of mesoscopic Josephson currents. It also follows that other 
quantities related to circulation as the helicity H = Jy u-wd?x may be non-zero 
too. The osmotic acceleration is given by 
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and the vorticity rate of change is 
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Defining the circulation as 
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and applying Stokes theorem we find 
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Hence, the circulation of the osmotic field does change with time and the Kelvin 
circulation theorem [10] [I7] is violated. It is important to note that as with the 
entropy production and the tensor equation (27), the non-zero vorticity effect 
is again a general result following from the inclusion of an imaginary quantum 
potential and does not rely on the phase assumption mentioned in section 2. It 
can be immediately verified that the same result holds by substituting the exact 
expression (10) in (30). 


4 Concluding remarks 


We have seen that the addition of an imaginary Bohm potential generates quan- 
tum friction and leads to certain interesting implications. It induces decay ef- 
fects in the amplitude of probability, modifies the exponential decay law, leads 
to entropy production and generates vorticity effects for the osmotic velocity. 
The phase of the wavefunction, and its derivative velocity field, has particular 
importance in the present formalism since it guides the osmotic velocity field. 
Furthermore, still in the hydrodynamic formulation a generalized quantum hy- 
drodynamic propagator has been obtained. In the domain of standard dynam- 
ics, by combining the two hydrodynamical equations a single non-homogeneous 
quantum Lighthill wave equation, manifestly time-symmetric, has been derived 
which in general admits both retarded and advanced solutions. The wider phys- 
ical implications of this result suggest that we may need to choose between 
backward causation and exploring the possibility of extending the standard for- 
malism beyond the unitary regime. 


Appendix 


Biparticle compound 


We consider a system of two spinless interacting particles in one dimension 
(e = 0). Following the derivation in [7] for the single particle case, we derive the 
probability density local continuity equation (p = p12(£1, £2), S = S12(£1, £2)) 
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and the corresponding Hamilton-Jacobi one 
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Taking the appropriate spatial derivatives we find 
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Using the formula 
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we obtain these Navier-Stokes equations for the momentum fields 
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As earlier we differentiate (38) over time and, (43) over x; and (44) over x2 and 
subtract from the first the sum of the two latter. Using the relation 
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The last terms in (46) or (47) indicate the non-separability of the biparticle 
compound. They exhibit particular conceptual interest as the operators act at 
two distant sites and cannot be ascribed to one or the other individual particle. 
The interparticle expressions for the many-body case follow this same pattern. 
Using index notation we write two-particle analogue of (26), setting € = 0, as 
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upper indices denote number of particle and lower the velocity field compo- 
Pp y- 
nents.) 
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